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The early thermalization and HBT puzzles in relativistic heavy-ion col- 
lisions studied at RHIC are shortly reviewed. The results of recent hy- 
drodynamic calculations that shed light on these two intriguing issues are 
presented. In particular, the role of the elliptic flow as a signature of early 
thermalization is critically examined. 

PACS numbers: 25.75.-q, 25.75.Ld, 24.10.Nz, 25.75.Nq 

1. Introduction 

The early thermalization and HBT puzzles in relativistic heavy-ion col- 
lisions were addresses simultaneously for the first time by Heinz and Kolb in 
2002 pQ (see also [2 [3] ) . In view of the first successes of relativistic hydro- 
dynamics applied to describe the hadronic spectra at RHIC [HESE], a short 
thermalization time scale of less than 1 fm (c = 1) as well as disagreements 
between the HBT data and the hydrodynamic predictions for correlation 
functions represented a serious challenge for theory. In the meantime, a sub- 
stantial progress in our understanding of those two puzzles has been made. 
At present, the agreement between the advanced hydrodynamic calculations 
and the HBT data is much better. However, the thermalization time scales 
assumed in such calculations remain short, well below 1 fm. This leaves still 
much place for a discussion of the early thermalization puzzle. The latter 
is most commonly solved by the assumption that matter formed in the rel- 
ativistic heavy-ion collisions at RHIC is a strongly interacting quark-gluon 
plasma (sQGP) 0[S]. 

As the notion of the "early thermalization puzzle" reflects simply our 
prejudice to accept an extremely short thermalization time (not supported 
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by the available microscopic calculations), the notion of the "HBT puzzle" 
is much more involved. First of all, the name "HBT analysis" is not very 
much informative. It refers mainly to the study of pion correlation func- 
tions (the two-particle distributions of identical pions with small relative 
momenta) , for which the names "pion interferometry" or "pion femtoscopy" 
seem to be more adequate. The term "HBT" emphasizes the analogy to the 
well known technique of second order intensity interferometry developed by 
Hanbury-Brown and Twiss [TO] to measure stellar angular sizes [Jj. The 
measurements of the pion correlation functions are quantified by the values 
of the "HBT radii" that define the range of the correlation functions. In 
this context, one should remember that the HBT radii i? s idci -Rout, and R\ ong 
give information about the space-time sizes of the regions where pions are 
correlated (so called homogeneity lengths [14] ) rather than about the actual 
sizes of the whole system. The persistent differences between the measured 
and modeled HBT radii were termed the "(RHIC) HBT puzzle". 

As was stated above, the two puzzles emerge in the context of rich ap- 
plications of relativistic hydrodynamics to describe heavy-ion collisions at 
RHIC. Contrary to the calculations based on hydrodynamics or kinetic the- 
ory, the one- and two-particle distributions may be consistently described in 
terms of simple, blast- wave-type models which assume appropriate physical 
conditions (temperature, transverse flow) on the properly chosen freeze- 
out hypersurface |15|. I16j. The problem remains, however, to validate such 
freeze-out conditions with the help of hydrodynamic or kinetic models which 
are regarded as more fundamental frameworks. 

More information and details about the HBT analysis and the HBT 
puzzle the reader may find in the recent reviews |17|. I18j . In this paper 
we discuss the HBT puzzle in the context of very recent hydrodynamic 
calculations. 



2. Resolving the HBT puzzle 

The common present opinion is that several improvements done in the 
hydrodynamic models may improve their predictions in such a way that the 
HBT puzzle is practically eliminated I2U] (discrepancy between the data 



This analogy is very often misunderstood. The particle momentum correlations are 
confused with the spacetime HBT correlations. Although the two types of the corre- 
lations have the common quantum statistical origin, the momentum correlations of 
identical particles yield the spacetime picture of the source, whereas the spacetime 
HBT correlations (dependence of the number of coincident two-photon counts on the 
distance between two detectors) provide the information on the characteristic relative 
three-momenta of emitted photons, which gives the angular size of a star without the 
knowledge of its radius and lifetime [111 1121 [T5] , 
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and theory remains at the level of about 10%). The discussed improvements 
include: 

i) the use of a realistic QCD equation of state |21U22j (instead of a simple 
bag equation of state), 

ii) early start of hydrodynamics (with the starting times as short as 
r eq ~ 0.10-0.25 fm [201 [23]), 

hi) inclusion of the pre-equilibrium flow [24, 25J, 

iv) inclusion of the shear viscosity [20] , 

v) accounting for the effects of fluctuations of the initial eccentricity [26] 
(the use of the participant plane rather than the reaction plane in the 
calculations of the elliptic flow), 

vi) the use of two-particle methods for the calculation of correlation func- 
tions [27] . 

vii) modified initial conditions (for example, the Gaussian profiles of the 
energy density in the transverse plane [19J), 

viii) fast freeze-out (for example, the use of a single-freeze-out scenario |28] 
where the collisions between the emitted hadrons are neglected), 

ix) inclusion of a complete set of hadronic resonances [27] . 

In our opinion the consensus has been reached for the points: i), ii), 
v), and vi). Moreover, the point ii) may be treated, to some extent, as 
equivalent to hi) and iv). The points vii)-ix) were proposed and emphasized 
inRef. ED. 




Let us now briefly discuss the physical arguments standing behind the 
points i)-ix): The use of the realistic QCD equation of state, e.g., see [2T1I22"] . 
which is much stiffer than the bag equation of state with large latent heat 
leads to less extended outward dimensions, which lowers the i? out /i? s ide ra- 
tio and makes it more compatible with the experimental data indicating 
R ov x/R s \de ~ 1- Similar effects are caused by an earlier start of hydrody- 
namics or by the inclusion of the pre-equilibrium flow [244 125] in standard 
initial conditions. For example, to avoid problems with choosing the appro- 
priate form of the initial transverse flow at r = 1 fm, one may decide to 
start the hydrodynamic evolution at an earlier time, r < 1 fm, and treat 
the results at r = 1 fm as new initial conditions with transverse flow. 

The presence of shear viscosity also increases the explosiveness of the 
collision |29t [5D| 131] . This can be seen by considering viscous corrections 
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to the stress-energy tensor. At the beginning of the evolution the velocity 
gradient is mainly longitudinal, which affects the stress-energy tensor by 
increasing the transverse pressure and decreasing the longitudinal pressure. 

It should be emphasized that the solution of the HBT puzzle cannot 
be reduced only to the correct reproduction of the HBT radii in a certain 
hydrodynamic calculation. In addition, one should describe well also the 
transverse-momentum spectra and the elliptic flow of hadrons. In other 
words, the acceptable solution of the HBT puzzle requires that both the one- 
and two-particle observables are uniformly reproduced @. In this context it 
is important to establish the most reliable method of the calculation of the 
elliptic flow in the hydrodynamic calculations. 

In Ref. [32] the important point was made concerning the inclusion of 
fluctuations of the eccentricity characterizing initial distributions of matter 
in the transverse plane. Such fluctuations increase the value of the elliptic 
flow. The inclusion of the fluctuations in the theoretical calculations helps 
to reconcile the predictions of the models with the data [32J. It was also 
shown by the same group [26] that the fluctuations improve the agreement 
between the HBT data and hydrodynamic calculations. 

The fluctuations of the eccentricity may be taken into account in a simple 
way by fixing centrality class and preparing "typical hydrodynamic initial 
conditions" for this class. This can be achieved, for example, by doing 
Monte-Carlo simulations of the Glauber model with GLISSANDO [33]. In 
this procedure the fluctuations are included in an averaged way. Definitely, 
a more basic but also a more time-consuming approach requires that the 
hydrodynamic evolution is applied after each generation of initial conditions 
by the appropriate model of the early stage (event-by-event hydrodynamics) . 
Such strategy is used now in Refs. [34"1 [35] . 

Early hydrodynamic calculations of the HBT radii used frequently sim- 
ple formulas relating -Rside, -Rout, and i?i on g to the space-time extension 
of the produced system. Nowadays, much more sophisticated methods are 
applied that try to mimic more closely the experimental situation. The two- 
particle methods are used in the framework of the event generators, e.g., 
THERMINATOR [36j. The correlation functions are constructed with the help 
of the emission function and the squared wave functions of the pion pairs 
|27j . If the wave functions include the Coulomb interactions, the Bowler- 
Sinyukov method [371 EE] is used to extract the HBT radii. The explicit 
calculations using simple parameterizations of freeze-out show that the ap- 
proach using the standard fitting without the Coulomb corrections agrees 
very well with the approach that includes the Coulomb interactions and 

2 The hydrodynamic predictions make sense in the soft hadronic sector where 
px < 1 — 2 GeV. Thus, the HBT puzzle discussed here refers to the soft hadronic 
observables. 
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uses the Bowler-Sinyukov method [271 [39]. 

In Ref. [19] a very much successful description of the soft hadronic 
observables studied at RHIC was achieved in the hydrodynamic approach 
using modified initial conditions for energy density in the transverse plane. 
Instead of using the Glauber model in the optical limit, the Monte-Carlo 
Glauber model [33] was used to determine typical initial conditions (distri- 
butions) for each centrality class. Then, the Gaussian fits were performed 
to those distributions, which were subsequently used as the input for the 
hydrodynamic calculations. It turns out that the initial Gaussian distri- 
butions lead to the freeze-out conditions which reproduce very well soft 
hadronic observables. 

At present no microscopic explanation of the Gaussian initial conditions 
is known. They do a good job because steeper profiles of the initial energy 
density in the transverse plane lead to larger pressure gradients at the center 
of the fireball and speed up development of the transverse flow. In this 
sense, the point vii) leads to similar effects as the points ii), iii), and iv). 
In addition, the initial Gaussian profiles lead to the Hubble flow [40] which 
is a kind of optimal flow for correct reproduction of the HBT radii in the 
blast- wave models [23] . 

The calculations pj|] used the concept of a single freeze-out [28]. In 
this approach the hadronic rescattering in the final state is neglected. To 
check if this assumption makes sense, the authors of []15] made estimates 
of the number of pion and proton collisions after freeze-out. This number 
was estimated as follows: pion straight-line trajectories after freeze-out were 
assumed and the number of collisions was counted, i.e., the encounters with 
other particles closer than the distance corresponding to the cross section. 
Averaging over all pions yields the number of these trajectory crossings 
about 1.5-1.7 per pion. These estimates indicate that the single-freeze-out 
scenario may be well accepted for pions. For protons, similar estimates are 
worse indicating non- negligible effects that may change the proton V2, see 
our discussion below in Sect. [5j 

In the calculations of the HBT radii it is important to include the effects 
of the resonance decays which make the studied systems larger by about 1 fm 
[27] . Most of the publicly available codes that simulate resonance decays 
use the resonance tables from SHARE [41] . 

3. Early thermalization problem 

Within the concept of a strongly coupled quark-gluon plasma short ther- 
malization times appear in the natural way, as the strongly interacting sys- 
tem equilibrates very fast. However, it is still debatable if the plasma in the 
early stage is indeed strongly coupled. In the last years many different ap- 
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proaches were used to discuss the problem of very early thermalization and 
very often such approaches considered the plasma as a weakly interacting 
system: 

i) The equilibration problem was studied within the parton cascade 
model. Initially, this model included only binary collisions but further 
developments took into account the gluon radiation in the initial and 
final states [12]. More recently, the advanced numerical codes have 
been developed, that stress the role of the multi-particle processes. In 
Ref. [13] the production and absorption 2 f> 3 processes are taken 
into account, whereas in Ref. [H] the three-particle collisions 3 o 3 
are studied. Within both approaches the equilibration is claimed to 
be significantly sped up when compared to the equilibration driven by 
the binary collisions. 

ii) One usually assumes that the initial partons are produced by hard 
or semi-hard interactions of partons in the incident nuclei. In the 
"bottom- up" thermalization scenario the initial state is described 
by the QCD saturation mechanism that is incorporated in the frame- 
work of the color glass condensate. Thus, the initial state is dominated 
by the small x gluons of transverse momentum of order Q s (one ex- 
pects Q s ~ 1 GeV for RHIC, and Q s ~ 2-3 GeV for the LHC 05]). 
The calculations performed within the "bottom-up" thermalization 
scenario [45] , where the binary and 2f}3 processes are taken into 
account, give an equilibration time of at least 2.6 fm |46] . see also 

iii) The equilibration of the system is most commonly understood as an 
effect of parton rescattering. An interesting phenomenon occurs, how- 
ever, that the equilibration is speeded-up by instabilities generated in 
an anisotropic quark-gluon plasma [48 [ 149] . This is so because the 
growth of the unstable modes is associated with the isotropization of 
the momentum distributions, which helps to achieve the full equili- 
bration. 

iv) It is also argued that the very process of particle production leads 
to the equilibrium state without any secondary interactions. For in- 
stance, Refs. |50[ [51] refer to the Schwinger mechanism of the particle 
production in the strong chromoelectric field. This approach explains, 
however, the equilibration of the transverse momentum only. The ap- 
proach of Refs. [S2J [S3], where the longitudinal momentum is also 
thermal, evokes the Hawking-Unruh effect. 
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v) An interesting physical scenario was also formulated, where the ther- 
malization is not an effect of collisions but a consequence of the chaotic 
dynamics of the non-Abelian classical color fields, coupled or not to 
the classical colored particles [54J. 

One may summarize different approaches to the early-thermalization prob- 
lem with the conclusion that newly developed simulations based on pertur- 
bative QCD as well as new non-perturbative frameworks may explain the 
thermalization times of about 1 fm. The problem remains to explain much 
shorter thermalization times used in the hydrodynamic codes that success- 
fully describe the data. For a more detailed discussion of these points see 

4. Hydrodynamic models with delayed thermalization 

The support for the idea of a very short thermalization/equilibration 
time at RHIC comes mainly from the observation of the large elliptic flow. 
This effect has been explained by the hydrodynamic (perfect-fluid) expan- 
sion with an early starting time. Recent hydrodynamic calculations chal- 
lenge this point of view. It turns out that it is possible to connect the hy- 
drodynamic perfect-fluid hydrodynamic evolution with the pre-equilibrium 
stage that is not completely thermalized. The results of such calculations 
are compatible with the data as we discuss below. 

4-1- Early free- streaming stage 

An approximation to the early-stage dynamics in relativistic heavy-ion 
collisions consisting of the free streaming (FS) of partons followed by a sud- 
den equilibration (SE) to a thermalized phase which subsequently undergoes 
a hydrodynamic evolution was proposed by Kolb, Sollfrank, and Heinz in 
[56j . In the meantime, this approximation has been frequently considered in 
the modeling of the early stages of the nuclear high-energy collisions in the 
context of equilibration. In particular, it has been considered in an investi- 
gation of the isotropization problem by Jas and Mrowczynski [57], as well 
as elaborated in the context of the early development of flow by Sinyukov, 
Gyulassy, Karpenko, and Nazarenko \58\ 159]. 

The early hydrodynamic approaches overlooked the fact that for non- 
central collisions, where the system develops spatial azimuthal anisotropy, 
an initial azimuthally asymmetric transverse flow may develop as a conse- 
quence of combined free streaming and sudden but delayed equilibration. 
Interestingly, the results of Ref. [40J obtained with and without the free- 
streaming stage are practically the same, see Figs. [TJand[5J We note that 
the initial conditions used in [40] are always Gaussian. 
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Fig. 1. The transverse-momentum spectra of pions, kaons and protons for the 
centrality class c=20-30% (upper panel) and the elliptic flow coefficient V2 for 
c=20-40% (lower panel). The darker (lighter) lines describe the model results ob- 
tained in Ref. [10] for the case with (without) free streaming. Data from [501 161) . 



In the case without free streaming, the hydrodynamic stage considered 
in [3D] starts at r = tq = 0.25 fm and the initial central temperature is 
Ti = 460 MeV for c=20-30% and Tj = 500 MeV for c=0-5%. In an alterna- 
tive scenario, the free-streaming starts at r = to = 0.25 fm and at r = 1 fm 
a transition to perfect-fluid hydrodynamics is made. The transition is de- 
scribed with the help of Landau matching conditions. In this case the central 
temperature at r = 1 fm is T, t = 305 MeV for c=20-30% and Tj = 330 MeV 
for c=0-5%. 

4-2. Early transverse hydrodynamics 

The energy-momentum tensor obtained from the free-streaming stage 
matches very smoothly to the form of the transverse hydrodynamics, where 
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Fig. 2. The pion HBT radii i? s idc , -Rout , -Riong, and the ratio i? ut/-Rside for central 
collisions. The darker (lighter) lines describe the results obtained in Ref. [10] with 
(without) FS+SE. The data from [52]. 



the longitudinal pressure vanishes |63t I64j . This behavior indicates that 
the free-streaming stage and the perfect-fluid hydrodynamics stage may be 
separated by a phase that is treated as an anisotropic fluid whose trans- 
verse pressure is much larger than the longitudinal pressure. Such pressure 
anisotropy appears also in the string models |50t [51~] where it is a conse- 
quence of the specific production mechanism. 

Having in mind very large pressure anisotropics present at the early 
stages of heavy- ion collisions, a model was constructed where the initial 
stage described by the transverse hydrodynamics is matched to the stan- 
dard perfect-fluid hydrodynamics [65]. In analogy to the model described 
in Sect. 14.14 the phase described by the transverse hydrodynamics lasts 
between r = 0.25 fm and r = 1 fm, and the transition to the perfect-fluid 
regime is described again by the Landau matching conditions. However, 
in contrast to the previous model the initial conditions are taken from the 
Glauber model in the optical approximation. 

The model results presented in Figs. [TH2] and [3H51 describe reasonably 
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Pt [GeV] 



Fig. 3. The comparison of the model [65] and experimental transverse-momentum 
spectra of pions, kaons, and protons. The data are taken from Ref. 60 . The 
model includes an early stage described by transverse hydrodynamics. The model 
parameters are defined in |65) . 




well the basic soft-hadronic observables measured at RHIC. The model with 
the transverse-hydrodynamics stage describes the data at the level of 20% 
(note the quantity A, displayed in Figs. EH5J which gives the relative dif- 
ference between the data and the model results). The model with the free- 
streaming stage describes the data even better (because of the Gaussian 
initial conditions used for hydrodynamics). 

Our two examples of the non-equilibrium early dynamics that precedes 
the perfect-fluid stage indicate that the early thermalization and the HBT 
puzzles may be circumvented in the pionic sector; the transverse-momentum 
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spectra, the elliptic flow coefficient x>2, and the HBT radii of pions are very 
well described in the models that do not assume early thermalization. In 
the two considered cases, the complete thermalization may be postponed to 
the comfortable time r = 1 fm. This value seems to be acceptable from the 
point of view of the recent microscopic calculations, see our comments in 
the end of Sect. 

Of the special interest is the fact that the elliptic flow of the pions is 
reproduced very well together with other observables. This result indicates 
that the large pion elliptic flow is not a good signature of very early ther- 
malization. 

The free-streaming stage as well as the transverse-hydrodynamics stage 
change into the phase treated as perfect fluid. In Refs. [BBJ |67] a physical 
scenario was considered where the transverse hydrodynamics was followed 
directly by a sudden isotropization and freeze-out (without the extended 
perfect-fluid stage). The results of Refs. [Ml EZ] are similar to those pre- 
sented in Figs. [3HHJ The approach with sudden isotropization and freeze-out 
gives slightly smaller HBT radii, on the other hand, their hp dependence is 
reproduced better, yielding a very good description of the ratio -R ut/-R s idc- 
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The similarities between the results presented in this paper and the results 
of Refs. [661 EZ] indicate that the details of modeling the equilibration 
transition may be not important. 

5. HBT or V2 puzzle? 

The largest discrepancies between the model predictions presented in 
Sects. 14.11 and 14.21 and the experimental data can be seen for the proton 
elliptic flow. The latter is probably too large because hadronic interactions 
in the final state and/or shear viscosity have been neglected. We know, 
however, that such interactions spoil the agreement with the HBT data. 
Clearly, some part of the "HBT puzzle" remains if we study pions, kaons, 
and protons. 

If we consider hydro dynamic evolution governed by the equations of 
perfect fluid and calculate physical observables at different freeze-out tem- 
peratures Ty, we find that with lowering of Ty the spectra become flatter, 
the splitting of the elliptic flow for different hadrons increases, and the HBT 
radii become too large. For large Ty (TV > 140 — 150 MeV, early freeze-out) 
the spectra are reproduced well, there is no V2 splitting, and the HBT radii 
are too small. On the other hand, at small TV (Ty < 140 — 150 MeV, late 
freeze-out) the spectra are too flat, the correct i>2 splitting is reproduced, 
and the HBT radii are too large. This type of behavior is the main source 
of the difficulties encountered in attempts to obtain a uniform description 
of soft-hadronic observables at RHIC. 

The optimal results are obtained for moderate temperatures, where the 
spectra are a bit too flat, a small V2 splitting is obtained, and the HBT 
radii are well reproduced. The points i)-ix) discussed in Sect. [2] improve 
significantly this type of optimization, with the effects presented in Sects. 
IPandlPl 

A real progress in improving the results of hydrodynamic calculations 
has been recently made by P. Bozek who included the effects of bulk viscos- 
ity. The inclusion of bulk viscosity lowers the pressure, makes the evolution 
longer, and leads to the V2 mass splitting. At the same time, the other 
hydrodynamic predictions describe well the spectra and the HBT radii [68J. 
Another recent calculation by P. Bozek shows that the early thermalization 
is required to reproduce well the directed flow [69J . This may suggest a shift 
of the ten-years-old paradigm: the directed flow rather than the elliptic flow 
is a signature of early thermalization. 

6. Conclusions 

In this paper we have tried to summarize recent attempts to solve the 
early thermalization and HBT puzzles at RHIC. Many improvements done 
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in modeling of the QCD equations of state, hydrodynamic initial conditions, 
as well as inclusion of dissipative effects improved the agreement between 
the hydrodynamic predictions and the data. Such improvements eliminate 
most of the discrepancies connected with the notion of the HBT puzzle. 
The combinations of inviscid hydrodynamics with non-equilibrium models 
of early stages suggest that the early thermalization puzzle connected with 
the elliptic flow may be circumvented. However, the new analysis of the 
directed flow [UU] shed completely new light on this problem. 
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